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ABSTRACT 

Aims. The observed variations in profiles of the interstellar aromatic infrared bands correlate with the object type and 
are indicative of PAH populations existing in different sources. Spectroscopic studies on PAHs can provide tools for the 
interpretation of variations accompanying the AIBs. As the observed spectra results from a mix of possible species in 
the region attempt is made to model this composite spectra by co-adding emissions from PAHs in different size groups. 
Methods. Theoretical IR data of PAHs having 10 to 96 carbon atoms is used to obtain emission spectra. The models 
are taken in size groups making up of small, medium and large PAHs. 

Results. The models show good profile match with observations for the 7.7 fim complex having sub-features at 7.6 and 
7.8 /xm. The 7.6 fim sub-feature dominates in the spectra of medium sized PAH cations matching observations from 
UV rich interstellar environments. The 7.8 /im component is more intense in the spectra of large PAH cations (model 
HI) correlating with observations from benign astrophysical regions. A possible interpretation for the observations of 
C — H out-of-plane bend modes and the weak outliers on the blue side of the intense 11.2 band is proposed. The 
models provide pointers to possible PAH populations in different regions. 



1. Introduction 

The mid-infrared emission bands at 3.3, 6.2, 7.7, 8.6, 11.2 
and 12.7 ^lm (3030, 1610, 1300, 1160, 890 a nd 79 cm^^) 
and b eyond (Gillett et a l. [19731 Cohen et al. [19891 Geballe 
et al. ITMSl ISO results [TM51 Peeters et al. [MHl Peeters 
et al. 1200 4a) observed in diverse astrophysical sources are 
popularly known as aromatic infrared bands (AIBs) . These 
bands have been attributed to emission from isolated gas 
phase PAH molecules present in the interstellar medium 
(ISM) (Leger fc Pu get [T9841 Allamandola et al. il985l 
Puget & Leger dnHSl Allamandola et al. nMS]) . The ubiq- 
uity of AIBs make PAHs an important interstellar family 
of molecules having significant carbon reserves. The emis- 
sion process of a PAH molecule involves the absorption 
of an energetic photon ranging from near UV to near IR 
(Allamandola et al. [T9891 Li & Draine [20021 Smit h et al. 
[MHl Mattioda et al. '2005a| Mattioda et al. l2005b|) . which 
renders it to a vibrationally excited state through internal 
vibrational redistribution (Allamandola et al. 119891 Puget 
& Leger [155^ . PAH molecules de-excite mainly through IR 
fluorescence with the emission bands corresponding to fun- 
damental modes of vibrations set up within the molecule. 
PAHs may also relax via visible phosphorescence and/or 
fluorescence and are proposed as carriers of the extended 
red emission (ERE) (Witt et al . [2006)) and the blue lumi- 
nescence (BL) (Vijh et aL [MI51) . 

Extrapolating the spacing between 6.2 and 7.7 /im 
PAH bands p roduced from experimental data, Hudgins & 
Allamandola |1999 bP have concluded that PAHs with 50- 
80 carbon atoms dominate the mid-IR emission. Studies 
testing the photo-physical stability of PAHs (Schutte et al. 
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[T993! AUain et al. I1996al AUain et al. I1996b| LePage et 
al. 2003) have reached a similar conclusion. Devoted ex- 
perimental (Szczepanski & Vala ll99^ Hudgins et al. 1994| 
Hudgins fc Alla mandola H MH H udgins & SandfordpSHl 
Cook et sl. mm Kim et al. lMm Kim & Saykally[2nni and 
theoretical research (Langhoff 1996; Bauschlicher fc Bakes 
2000; Bauschhcher [^0 2| Pat hak 2006, Pathak fc Ras togi 
2005t Pathak fc Rastogi [20061 Pathak & Rastogi [20071) has 
resulted in considerable IR information of large number 
of PAHs of varying sizes and ionized states. These results 
serve as inputs to theoretical emission models (Bakes et 
al. I2001al Bakes et al iSOOlbl Pech et al. J002; Jobhn et 
al. "2002} Mulas et al. l2003l Mulas et al. [M)6a Mulas et 
al. I2006bp for direct comparison with observations. Such 
modeling studies have emphasized the contribution of PAH 
cations towards AIBs along with neutral PAHs contributing 
to a lesser extent. The fractional abundance of cations and 
neutrals depends on the astrophysical environment where 
the AIBs are being observed. Current understanding points 
towards dominance of cations in harsh UV rich conditions 
of reflection nebulae and HII regions while a mixture of 
neutrals and cations in relatively cool environments of late 
type carbon rich stars (Allamandola et al. I1999P . 

Over the years, developments in observational accuracy 
using space telescopes have improved the quality of data 
and a large and diverse sample of astrophysical sources have 
been examined that exhibit these mid-IR bands (Peeters et 
al. [20021 Peeters et al. i^OlMal V an Diedenhoven et al.[2004|) . 
Observations from ISO (fT996l) and SPITZER (|2004l) have 
provided good quality data that has revealed significant 
variations and finer profile details of the individual AIBs. 
Of great interest is the peak position of 6.2 ixm band, the 
broad composite 7.7 band (Peeters et al. l2002p and vari- 
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ations accompanying the asymmetric 11.2 /im band, fea- 
tures at 12.7, 13.3 fim and at higher wavelengths (Hudgins 
& AUamandola 1999a; Hony et al. 2001). Few new features 
detected around 6.7, 10.1, 15.8, 17.4 and 19.0 fim have fur- 
ther added to the AIB numbers (Werner et al. |200 4: Peelers 
et al. I2004bp . These revelations make studies on different 
PAH groups important to obtain better understanding of 
AIBs. 

In this report we attempt to model the variations asso- 
ciated with AIBs using our existing IR database of PAHs 
(Path ak [20061 Pathak fc Rast ogi 2005, Pathak & Rastogi 
[2006 1 Pathak & Rastogi [20071) . The IR data consisting of 
PAHs with 10 to 96 carbon atoms is computed using the 
GAMESS ab-initio program (Schmidt et al. [I993|) under 
DFT using B3LYP functionals in conjunction with 4-31G 
basis expansion. To model the observed profile variations in 
AIBs first the thermal emission model (Schutte et al. 119931 
Pech et al. I2002p is used to obtain the emission spectra of 
each PAH. Then for PAHs in different size groups the emis- 
sion spectra is co-added to model AIB spectra. The compar- 
isons with observations give insight into possible PAH size 
distribution and their processing in different astrophysical 
regions. Attempt is made to put constraints on the PAH 
population residing in different interstellar environments. 

2. Emission model 

For a meaningful comparison with observations emission 
spectra of PAHs is needed (Schutte et al. ,1993; Cook & 
Saykally 1998; Pech et al. [200 1 Mu las et al. [2003|) . Leger 
et al. (1989) and Schutte et al. (|1993p have shown that cool- 
ing of a PAH via transitions in an emission cascade can be 
described by the thermal approximation if the average en- 
ergy of the mode under consideration is small compared to 
the total energy of the excited PAH. The thermal model is 
known to break down for emission in low energy modes near 
the end of the emission cascades ( Joblin et al. 120021 Mulas 
et al. [2I)06a; Mulas et al. I2006b[) . In the present work a 
strong background radiation field, corresponding to black- 
body temperature T — 40,000isr, is considered. This en- 
sures that the total internal energy of excited PAH is much 
higher than the energy of any emission mode and thermal 
approximation holds. A similar approach is taken by Cook 
& Saykally (1998) and Pech et al. f^UU^ . 

We compute the emission spectrum for each individual 
PAH with calculated absorption spectrum used as input. 
As in Schutte et al. (1993), the PAHs are considered in a 
exciting radiation field similar to that surrounding Ori 
C, corresponding to blackbody temperature T* = 40, OOOA'. 
The PAH molecules absorb photons of all energies with cut 
off at 13.6 eV giving rate of photon absorption as: 

(■13.6 TD 

Rabs = / (1) 

where Bi, is the Planck function, ai, is the photo- 
absorption cross-section and v are the absorbed UV ra- 
diation frequencies. For each PAH tr^ are taken from 
http://astrochemistry.ca.astro.it/database/ (Malloci et al. 
2007p . The differences in the shape of the photo-absorption 
is apparent mostly in the low energy part of the spectrum, 
which for the UV-rich spectrum in the present emission 
models is irrelevant. So, for PAHs having no entry in the 
database, scaled photo-absorption cross-section of PAHs 



closest in size (number of carbon atoms) is used. The scal- 
ing is approximated by a factor Nc2/N^ci (Mulas et al. 
2006b), with PAH having Nci carbon atoms as proxy for 
larger molecule with Nc2 carbon atoms. 

The absorption of a high energy photon of frequency v 
excites the PAH molecule to an internal energy equivalent 
to peak temperature Tp that depends on the heat capac- 
ity of the PAH. In the harmonic approximation this Tp is 
obtained from: 

^^^^ " S exp(/ic uj,/kT) - 1 

where i corresponds to vibrational modes within individ- 
ual PAHs having frequency iVi in cm~^ and m is the total 
number of vibrational modes (3N - 6; N being the number 
of atoms) . 

To ensure the validity of the thermal approximation fall 
in peak temperature for each PAH is calculated assum- 
ing integrated emission of highest frequency modes corre- 
sponding to the C — H stretch vibrations. For naphthalene, 
the smallest PAH in the sample, peak temperature falls by 
11.3%, which is maximum. While for CqqH24,, the largest 
PAH, fall in peak temperature is less than 0.5%. This is 
sufficiently small to justify the use of thermal model in the 
exciting radiation field considered. 

The molecule cools down from the peak temperature Tp 
through its different vibrational modes with cascade tran- 
sitions from levels w — > w — 1. The emission photon flux 
of the i*^ mode is given as: 

= a\° X [exp(/ic uj.,lkT) - Vf^ (3) 

where the Einstein coefficients A^'^ are obtained from the 
absorption intensities Si (in units of Km/mol) (Cook & 
Saykally dHMI) by the relation: 

A\'^ ^ (1.2512 X 10"^) Luj 5, (4) 

For a fall in internal energy by AC/ the fractional energy 
emitted in the i*'' mode is then given as: 

A£;,(r) = ^tl"' X Ac/(r) (5) 

The fractional energy Ei is integrated over the cooling range 
from Tp to a temperature of 50 K below which the energy 
emitted is negligible. At/ corresponding to fall in temper- 
ature by 1 K is taken at each T. This emitted energy is 
weighted by the rate of photon absorption (Eq. 1) and in- 
tegrated over the whole distribution of absorbed photons. 

The emission spectra thus obtained for each PAH is then 
used to model composite spectrum from groups of PAHs. 
Spectra of PAHs in each group are co-added assuming equal 
number of each specie in the ISM. 

3. Spectral models and their astrophysical 
significance 

AIBs result from composite emission of a number of PAHs 
consequently profile variations associated with AIBs in dif- 
ferent sources reflect the presence of different PAH groups 
therein. Factors like anharmonicity and rotational broad- 
ening also affect line widths and shape (Pech et al. 2002]). 
Anharmonicity induces small asymmetry in the profile and 



A. Pathak & S. Rastogi: Modeling AIBs with co-added spectra of PAHs 



3 



rotational widths too are small and get smaller for larger 
PAHs. Contribution of these effects is < 5 cm~^ (Mulas 
et al. [2006 a: Pech et al. 2002), while typical peak width 
under interstellar conditions is 30 cm~^ (AUamandola 
et al. I1989p . Experimental measurements on a few individ- 
ual PAHs show a near linear dependence on band width 
and temperature with intrinsic widths between 10 - 20 
cm-i around 1000 K (Joblin et al. lTMSl Pech et a\. ?KM\i . 
In obtaining the composite emission spectra we neglect 
anharmonicity and rotational broadening and assume a 
Lorentzian profile with FWHM of 20 cm~^ . Several modes 
closer than this may collocate together and result in a single 
broad peak. Three distinct models are presented to under- 
stand the gross changes (intensity as well as profile varia- 
tions) in the PAH spectra with different PAH size groups. 
Individual PAHs used in the models are displayed in Fig. 

m 

Model I comprises of co-added emission spectra of small 
PAHs (< 20 carbon atoms) consisting of naphthalene, an- 
thracene, phenanthrene, pyrene, tetracene, chrysene, 3,4- 
benzophenanthrene, triphenylene, perylene and benzopy- 
rene. Model II (PAHs with 20 to 40 carbon atoms) is derived 
from a mixture of anthanthrene, coronene, bisanthene, ova- 
lene, C28H14, C30H14 and CagiJie- Model HI uses the co- 
added spectra of large PAHs (> 40 carbon atoms) which in- 
cludes C4sHis, C^tHiq, CQ2H20, CqqH20, CS0H22, CgoH24 

and CgeH24. All the three models consider equal number of 
constituent PAHs. Emission spectra of the three models as 
neutral and cation are compared in Fig.[2]and Fig. [31 In the 
composite spectra each peak may result from the overlap 
of different modes in different PAHs. Therefore only a gen- 
eral assignment of intensity peaks with vibrational modes 
is possible. 

The spectra of neutrals (Fig. ^ is dominated by the 
C — H stretch and C — H out-of-planc bend vibrations. 
Comparison of the absorption calculations with experimen- 
tal data wherever available shows that the C — H stretch 
intensity is overestimated by about 1.4 to 2.0 times (Pathak 
& Rastogi [2007) . Use of larger basis sets, i.e., 6-31G or 6- 
31G**, in place of 4-31G, provide better representation of 
the orbitals and consequently a good C — H stretch inten- 
sity match with experiments. Use of larger basis for some 
of the PAHs incorporated in model I shows that computa- 
tional effort increases greatly and the scaling procedure also 
gets complicated but there is negligible change in intensity 
or position of modes other than C — H stretch (Pathak & 
Rastogi [^007p . Therefore, for the study of profiles and rela- 
tive intensity of all other modes the 4-31G basis is suitable 
while giving qualitative information of the C — H stretch. 

In the cascade emission model low temperatures con- 
tribute little towards the intensity of modes at higher fre- 
quencies. The contribution is diminished further in case of 
large PAHs that are excited to lower peak temperatures. In 
large PAH model III the intensity oi C — H stretch mode 
is small compared to the models with smaller PAHs. The 
observed C — H stretch (3.3 /zm) and C — H out-of-plane 
bend (11.2 jim) intensity ratio is a suitable parameter con- 
straining the size distribution of PAHs (Schutte et al. 119931 
Pech et al. [M7^ . 

The absorption intensity oi C — H stretch vibrations re- 
duce drastically upon ionization while small intensity vari- 
ations are observed for the C — H wag modes. Intensity of 
C — H stretch mode depends on the charge density near the 
hydrogen atoms (Hudgins et al. 120011 Pathak & Rastogi 



120051 Pathak & Rastogi I2006P . In cations this (positive) 
charge density increases consequently the C — H stretch in- 
tensity gets drastically reduced. In larger PAH cations due 
to the distribution of the acquired positive charge over more 
number of atoms the absorption in C — if stretch mode 
is significant (Pathak & Rastogi 120061 Pathak & Rastogi 
I2007P . This and the temperature effect in cascade emission 
model is reflected in the composite emission spectra (Fig. ^ 
showing significant C — H intensity in model II and reduced 
intensity in model I and model HI. 

3.1. C — C stretch modes 

Among the AIBs the 6.2 and 7.7 fim emission features are 
most intense. The C — C stretch vibrations set up in ionized 
PAH molecules give rise to these bands. The 6.2 /im band 
has some contributions from C — H in plane bend modes as 
well (Pathak & Rastogi I2005P . Profile variations observed 
in these bands are a direct measure of the background en- 
vironments that excite the PAHs. 

Peelers et al. (2002) report detailed ISO observations 
of 57 diverse sources in the 6 to 9 jim range and based on 
intricate position and profile variations classify these bands 
into three classes A, B and C. Class A has the 6.2 /xm 
band maximum in between 6.20 and 6.23 /im and the 7.7 
/iTO complex peak at 7.6 ^m. Bands between 6.23 and 6.29 
jjLm and a peak at 7.8 fim are classified as class B. Bands 
centered at 6.29 and redder are classified as class C . Out of 
the 57 sources reported, 42 belong to class A, 12 to class B 
and 2 sources are classified as C . Similar profiles have also 
been reported by other observations for reflection nebulae 
and Herbig Ae/Be stars (Bregman & Temi [^0051 Sloan et 
al. 1^005) . 

All three PAH cation models (Fig. [3]) show intense broad 
peaks corresponding to 7.7 and 6.2 /im modes. The peaks 
are at 1330 (7.52 ^lm) and 1526 cm"! (6.55 /xm), 1347 (7.42 
^m) and 1537 cm~^ (6.50 ^m) and at 1285 (7.78 jim) and 
1564 cm~^ (6.39 urn) respectively for models I, II and III. 
The C — H in-plane vibrations have peaks around 1215 
cm~^ (8.23 ^m) in all the models with little variation in 
the peak position. 

Observations show that the 7.7 /im AIB incorporates 
two sub-components, at 7.6 (1315 cm~^) and at 7.8 
(1282 cm ~^) (C ohen et al. il989] Bregman et al. [TO5gi 
Peelers et al. I2002p . The 7.6 /xm feature tends to domi- 
nate the emission spectra of regions involved in process- 
ing of PAHs and other molecules i.e., star forming regions, 
HII regions and reflection nebulae, etc. The 7.8 /im com- 
ponent dominates in environments where the PAHs are rel- 
atively fresh and unprocessed for example, planetary neb- 
ulae and objects evolving from Asymptotic Giant Branch 
(AGB) phase of stellar evolution. Similar sub-components 
for the 7.7 /im band are obtained in the present models. 
This provides useful insight regarding the PAH size and 
charge state in different interstellar environments. 

Fig. [H shows the expanded spectra, in the 1250 - 1400 
cm~^ region, for the cation models of Fig. [3] The two par- 
allel lines at 1285 and 1315 cm~^ enclose the observed 7.7 
/im complex. Two sub- features which constitute the 7.7 
/im band having a separation of 30 - 40 cm~^ are evi- 
dent in the three models. Model I has a dominant feature 
at 1330 cm^^ (7.52 /xm) corresponding to the observed 7.6 
/im component along with a shoulder at 1300 cm~^ (7.7 
//m). In model II the lower wavelength component around 
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1347 cm^^ (7.42 fim), with shoulder at 1325 cm~^ (7.55 
/zm), dominates while the higher wavelength sub- feature 
at 1285 cm~^ (7.78 /im) is less intense. Profiles similar to 
model I and II, with stronger lower wavelength component, 
have been observed in UV-rich environments of HII regions 
and reflection nebulae and have been classified as A' pro- 
files (Fig. 13 in Peelers et al. ' 2002p . Medium sized compact 
PAHs with 40 - 50 carbon atoms, as present in model II, 
are quite stable and may survive in strong UV environments 
providing good representation of PAHs populating such in- 
terstellar sources. Observations also indicate that sources 
dominated by the 7.6 /im sub-feature always exhibit a dis- 
tinct 7.8 fj,m sub-feature (Peelers et al. [2002) . Model II 
fulfills this observational norm well. 

In model III that comprises large PAHs, the 7.7 
lim band is dominated by the higher wavelength sub- 
component at 1285 cm~^ (7.78 fim) with an apparent shoul- 
der around 1315 cm~^ (7.60 fim). Unlike model II it is not 
distinctly resolved but appears as a shoulder to the main 
peak. Observations show that sources with dominant 7.8 
/zm component do not always have a very clear and dis- 
tinctive 7.6 fj,m component (Peelers et al. [2D02; Bregman 
& Temi 120(15)) ■ The profile of model III conforms with these 
observations of relatively benign astrophysical regions clas- 
sified as B' profiles (Peelers et al. I2002p . 

For direct comparison and matching with observations, 
intensity ratios of the two components of the 7.7 yCtm com- 
plex and their wavenumber separation in the three models 
is presented in Tabic [TJ Similar ratios from observations 
by Peelers et al. (j2002p have also been provided. A good 
match in terms of intensity ratio is present with different 
observed sources. The profile observed in model I is simi- 
lar to the classification A' . Model II shows a separation of 
almost 62 cm~^ between the two components of the com- 
posite 7.7 band with an intensity ratio of 1.06. The 
separation is a bit more than that is observed for A' pro- 
files. The intensity ratio is also less than what is measured 
for NGC 2023 but is closer to that in IRAS 23133. Similar 
match has been shown to exist with the co-added spec- 
tra of medium sized pericondensed PAH cations (Pathak & 
Rastogi [20061 . Model HI shows a clear agreement with the 
observed B' profiles. The separation between the two com- 
ponents of the 7.7 band as well as the intensity ratio 
match well especially for NGC 7027. 

The region of 6.2 fxm feature, i.e. 1450 to 1650 cm~^ 
region of Fig. [31 is expanded in Fig. [5l There is no position 
match with the 6.2 (1610 cm^^) AIB in any of the three 
models. The co-added spectra of small PAHs (model I) has 
this band centered at 1526 cm~^ (6.55 /xm) which is far off 
from the observed 6.2 fj,m band. As the PAH size increases a 
blue shift, similar to available experimental data (Hudgins 
& AUamandola 1999b)) . is observed. The co-added spectra 
of medium sized PAHs (model II) has this feature at 1537 
cm~^ (6.51 /im) and for large PAHs (model III) only a small 
additional shift to 1564 cm~^ (6.39 /im) is present. Few 
PAH cations eg., C28-ffi4, Cas-ffie and CgoH24 (Pathak & 
Rastogi [20061 Pathak & Rastogi [2007)) independently show 
intense feature around 6.29 (1590 cm~^) corresponding 
to the observed class C band (Peelers et al. , [2002() . These 
PAHs contribute to the small shoulder observed beyond 
1590 cm-i (6.29 fim). 

Failure to explain the position of the 6.2 fim AIB with 
co-added model spectra of pure PAHs points towards other 
possibilities. Recent theoretical calculations on nitrogen 



substituted large PAH cations (PANHs) proposes to ex- 
plain the variations in the peak position of the 6.2 fj,m in- 
terstellar emission feature (Hudgins & AUamandola 20031 
Hudgins et al. l2005p . Incorporating nitrogen inside the ring, 
i.e., replacing one or more carbon atoms results in a shift of 
the band from 6.4 fim to 6.2 fim. Similar blue shift is also 
reported for Iron - PAH coordination complex (Simon & 
Joblin [2l)07p . Hydrogenated, de-hydrogenated and multiply 
charged PAHs, PAH clusters, anions and substituted PAHs 
(Bakes et al. I2001al Halasinski et al. [20051 Malloci et al. 
[20051 Malloci et al. [20071 Vuong & Foing [2000]) also need 
careful attention and they cannot be completely neglected. 
Either these or a combination of these with pure PAH ions 
may be responsible for the 6.2 fim emission feature. 

3.2. C — H out-of-plane bend modes 

The 11.2 fim (893 cm~^) AIB represents solo out-of-plane 
wag modes in PAHs. The weak 12.7 fim (787 cm-i) band 
and features at longer wavelengths represent trio and quar- 
tet hydrogen wag modes. ISO observations by Hony et al. 
(|200ip conclude that the spectra of planetary nebulae and 
evolved carbon-rich stars, where PAHs are supposed to be 
synthesized and are still unprocessed, have a strong 11.2 fim 
feature, weak 12.7 fim band and features at longer wave- 
lengths. The situation is almost reversed in case of UV-rich 
environments where the 12.7 fim band is as intense as the 
11.2 fim band. 

The spectra of small (model I) and medium sized (model 
II) neutral PAHs (Fig.O have strong quartet hydrogen wag 
modes peaking around 750 cm~^. Mixing of duo and trio 
hydrogen wags result in a broad peak around 815 cm^^ in 
small PAHs. For PAHs comprising 20 to 40 carbon atoms 
(model II) the intense solo wag mode near 900 cm^^ (11.11 
fim) shows best position and intensity agreement with the 
observed 11.2 fim AIB. The spectrum of model HI has this 
feature blue-shifted to around 915 cm~^ (10.93 fim). Model 
HI also displays discrete less intense features due to duo and 
trio hydrogen out-of-plane bend vibrations at 840 and 800 
cm^^ respectively. 

Upon ionization the C — H out-of-plane bend modes in 
all the PAHs show significant shift in band positions but 
intensity variations are small (Fig.[0]). Model I has intense 
peak at 760 cm^^ (13.16 fim) and less intense feature at 864 
cm~^ (11.57 fim). These features in model I correlate with 
observations of weak bands present on the long wavelength 
side of the 11.2 /xm band (Witteborn et al. [TM9p . The co- 
added emission spectra of medium sized (model II) PAH 
cations has an intense feature at 920 cm~^ and significant 
peaks at 755 and 785 cto~^. Large PAH cations (model HI) 
have prominent solo hydrogen (C — H) out-of-plane bend 
modes producing an intense peak at 933 cm^^ and weak 
features between 800 and 850 cm~^. The shifts in band 
positions of hydrogen out-of-plane modes with ionization 
and increasing PAH size are apparent in the three models. 

Sloan et al. P^99p observed complex less intense struc- 
tures on the short wavelength side of the principal 11.2 
fim AIB. They referred them as "blue outliers" . These blue 
outliers were more noticeable near the central exciting ob- 
ject. The 11.2 fim feature adopts a smoother profile in ob- 
servations away from the central source with weak or no 
blue outliers. These observations correlate with our mod- 
els. Weak features on the blue side of the 11.2 fim have 
been pointed by arrows in Fig. [51 Existence of stable com- 
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pact medium sized PAHs or cations in the vicinity of the 
central source seem to be a proposition. This also correlates 
with the observation that UV rich environments where 6.2 
and 7.7 /im show A type profile also have A type 11.2 fim 
(Van Diedenhoven et al. l2004p . 

3.3. Intensity Ratios 

The intensity of the AIBs may not be directly correlated 
with spectral models based on theoretical or experimental 
measurements but the ratios of different features prove to 
be highly useful. Systematic variations associated with in- 
tensity ratios for different modes may be used to constrain 
PAH size, structure and ionization state in the ISM. In neu- 
tral PAHs, where the spectra is dominated hy C~H stretch 
and wag modes with weak C — C stretch mode intensities, 
the 6.2 and 7.7 ^.m to 11.2 fim band ratio is small but for 
cations tremendous increase in the value is present, which is 
an indicator of the ionization state of PAHs. The emission 
model intensity ratios of few features are presented in Table 
[2l Since the models fail to reproduce the 6.2 jim feature the 
presented ratios are not directly relevant for observations 
but point towards a trend. 

For cations, the value of 6.2 to 7.7 /xm band ratio is 
smallest for large PAHs. The ratio of 11.2 to 12.7 ^im PAH 
bands for neutral PAHs sheds light on the shape (peripheral 
structure) of PAHs (Hudgins & AUamandola I1999al Hony 
et al. l200ip . As the PAH size increases, the ratio appears to 
increase steadily. Since, small PAHs inevitably have more 
corners, leading to more duo and trio hydrogens, the inten- 
sity of their wag is strong. Compact and stable large PAHs 
considered in the models have straight edges resulting in 
increased number of solo hydrogens and the spectra has an 
intense 11.2 /im band. 

3.4. PAH formation and processing 

The correlation of observations with model spectra, from 
collective emission of PAHs, presented in this report pro- 
vides leads to understanding the formation and processing 
of PAHs in the ISM. The correlation of 7.7 B' profile 
with model HI suggests formation of large PAHs (having 
around 100 carbon atoms) in outflows of late-type carbon 
rich stars. Since, ionization energy for large PAHs is less, 
they are more likely to be cations. The spectra of these 
large PAH cations dominate benign regions of the ISM. 
Observations of the 6.3 fim (1585 cm"^) AIB (C class; 
Peeters et al.. l2002p in these regions suggests that the PAHs 
initially formed are pure and unsubstituted. The process- 
ing of these large systems starts as soon as they are formed. 
Ruled by the ISM conditions equilibrium is established be- 
tween formation and destruction. A likely scenario is that in 
harsh ISM conditions breaking and reduction of large struc- 
tures along with multiple ionization of medium sized PAHs 
dominates. Observation of A' profiles in UV-rich environ- 
ments of ISM suggests the abundance of compact medium 
sized PAH cations. This is also confirmed by the observa- 
tion of a strong 12.7 /im band in UV-rich environments. 
Formation of multiply charged PAHs, substituted PAHs 
(PANHs), PAH complexes and clusters is suggested by the 
observation of 6.2 /tm A class profile in harsh strong UV 
flux environments (Peeters et al. [^02p . 



Taking cue from these models more specific modeling 
for particular observations can be made by varying the 
concentrations of the individual PAHs. These can provide 
better understanding of the chemistry of the region. A 
combination of PAHs from model I and II with differ- 
ent concentrations is taken to compare with observed 
7.7 /xm feature of reflection nebula NGC 2023. The 
most suitable combination obtained has Perylene(14%), 
Triphenylene(ll%), Benzo-pyrene(9%), Ovalene(25%), 
Bisanthene(9%), C3oi?i4(17%), C38i?i6(15%) and is shown 
in Fig. Elja). The lower wavelength component dominates 
the 7.7 /fm band while the higher wavelength sub-feature 
is slightly less intense. These have been classified as 
A' profiles (Peeters et al. I2002p observed in UV-rich 
environments of HII regions, reflection nebulae and other 
star forming regions. Starburst galaxies, where there is 
copious star formation activity, also have similar profile 
(e.g. average spectra of 11 Starbursts Fig.l Lutz et al. 
1998). The 7.7 /im band in galaxies is useful in the study 
of AGN-Starburst connection (Schweitzer et al. 120061 Tran 
et al. [200n Lutz et al. [T998)) . 

Fig.[7^b) compares the 7.7 jim feature of model III with 
post-AGB star HD 44179 that is classified as B' profile 
(Peeters et al. I2002p of relatively benign astrophysical re- 
gions where PAHs are supposed to be synthesized. A slight 
shift in the band position and an additional feature near 
1215 cm~^ (8.2 /im) mark the difference in agreement with 
observation. 

It is interesting to compare the models with observation 
over the whole range from 5-15 /im, as shown in Fig. [5) 
Though there is gross similarity, in both cases it is seen that 
apart from the 7.7 fim band other prominent features are 
either shifted or have intensity mismatch. The 6.2 /im band 
lies at a higher wavelength while the 11.2 /im band is blue 
shifted in the model spectra. Considering PAH population 
in ISM to be even partially similar to the present models, 
observations should show the strong features of the models 
e.g. around 6.5 and 10.8 /im. 



4. Conclusions 

AIBs are ubiquitously observed in diverse astrophysical 
sources with their origin commonly ascribed to interstel- 
lar PAHs. Recent space based observations have enabled 
high quality infrared data which has introduced variations 
in profiles of the observed mid-IR spectra. These variations 
correlate with the object type and are indicative of PAH 
populations existing in different sources. To gain insight 
into the profile variations of PAH bands and for their better 
interpretation three distinct models using PAHs of different 
size groups are presented. Correlation of these models with 
observations throws light on the formation and processing 
of PAHs in different astrophysical environments. 

The models show good profile match with observations 
for the 7.7 /im band. This feature is not a single peak but 
a combination of two sub-components at 7.6 and 7.8 /im 
that are also distinct in the presented models. The 7.6 /im 
sub-feature is found to dominate in the spectra of small 
to medium sized PAH cations matching with the observa- 
tions of UV rich interstellar environments (Fig.[7^a)). The 
7.8 /im component is more intense in the spectra of large 
PAH cations correlating with the observations of benign 
astrophysical regions (Fig. [7Kb)). This clearly indicates the 
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formation of large PAHs in the outflows of carbon rich stars 
that transform to medium sized ones upon processing. 

The models also correlate well with the observations of 
the C — H out-of-plane bend modes and the weak outliers 
on the blue side of the intense 11.2 /im band. Variation in 
the observed intensity of 12.7 and 11.2 bands points 
towards large PAHs in planetary nebula and medium and 
small PAHs in UV rich environments of HH regions and 
reflection nebula. 

Varying concentrations of the individual PAHs and con- 
sidering the strength of exciting radiation may allow mod- 
eling of specific astrophysical objects. But the feature corre- 
sponding to 6.2 /xm AIB is not correlated with the present 
set of PAHs. Also strong features of the models, around 
6.5, 8.2 and 10.8 /xm, should show up in observations if the 
PAH population in ISM is similar to those considered here. 
The models presented do provide pointers to possible PAH 
populations in different regions, yet for a complete and re- 
alistic model same group of PAHs should provide match for 
aU the AIBs. 

Failure to explain the complete profile of AIBs may be 
due to the incompleteness of the sample. A bigger sample 
of PAHs that includes substituted PAHs, PAH coordina- 
tion complexes, hydrogenated and de-hydrogenated PAHs 
and different charge and ionization states of all these may 
be required. Since laboratory data for most such systems 
may be difficult to obtain quantum chemical approach is 
most suitable, which needs to be refined by using larger 
basis sets. A better modeling approach should incorporate 
anharmonicity and hot band shifts. Some work has been 
done in this direction (Joblin et al. 11995! Pech et al. 120021 
Mulas et al. I2006al Mulas et al. :2006b; Cane et al. [M?7|) 
but exact identification of individual PAHs based on these 
is still difficult. 
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Table 1. Band position and strength ratios of the two com- 
ponents of 7.7 ^.m composite. 



Model 



7.6 



7.8 



r>m{cm-y Diff^fc, (em-^) Jt.s/Jt 



Model I 
Model II 



1330 
1347 



1300 
1285 



Model III 1315 1285 



30 
62 



30 



-30 



-45 



1.32 
1.06 



0.65 



^Difference of the band position and intensity ratio of the two 
components of the 7.7 fj.m composite as calculated from the 
theoretical spectral models.; 

■^The observed results are from Peeters et al. 120021 Table 2. 



Table 2. Intensity ratios" for the main PAH features. 



Model 


6.2/7.7 

(PAH+) 


6.2/11.2 

(PAH+) 


7.7/11.2 
(PAH+) 


11.2/12.7 
(Neutral PAH) 


Model I 


1.15 


7.65 


6.66 


0.44 


Model II 


1.22 


2.12 


1.74 


2.00 


Model III 


1.01 


1.70 


1.68 


4.39 



a -Ratios are for emission models as in Fig. [2] and Fig. |3] 




Fig. 1. PAHs constituting the three models. 
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Fig. 2. Model spectra of co-added emissions from neutral Fig. 3. Model spectra of co-added emissions from PAH 
PAHs with (I) less than 20 C atoms, (H) 20 to 40 C atoms, cations with (I) less than 20 C atoms, (H) 20 to 40 C atoms, 
(HI) more than 40 C atoms. (HI) more than 40 C atoms. 




Fig. 4. Emission spectra of three cation models in the 1250 
to 1400 cmr^ region. The observed 7.7 /xm AIB hes within 
the two paraUel hnes. 



Fig. 5. Emission spectra of three cation models in the 1450 
to 1650 cm^^ region. The observed 6.2 /xm AIB lies within 
the two parallel lines. 
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Fig. 8. Same as Fig. 7 in the range 5-15 ^m. (a) PAHs 
with less than 40 C atoms in varying proportions (see text); 
(b) more than 40 C atoms (model III). 



frequency {cm"') 



frequency {cm ') 



Fig. 6. Emission spectra of three models in the 600 to 1100 

cm^^ region. Arrows point to blue outliers of 11.2 /xm band. 
Dotted vertical lines represent central peak of 11.2 and 12.7 
Hm AIB. 



a b 




1400 1300 1200 1400 1300 1200 

frequency (cm ') frequency (cm"') 



Fig. 7. PAH cations emission model in the 1200 to 1400 
cm~^ region compared with observed spectra (a) PAHs 
with less than 40 C atoms in varying proportions (see text); 
(b) more than 40 C atoms (model HI). The 7.7 /zm band 
lies within the two vertical lines. 



